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A  COMPARISON  OF  THE  EXPERIMENTAL  SUBSONIC  PRESSURE  DISTRIBUTIONS  ABOUT 
SEVERAL  BODIES  OF  REVOLUTION  WITH  PRESSURE  DISTRIBUTIONS 
COMPUTED  BY  MEANS  OF  THE  LINEARIZED  THEORY1 

By  Clahence  W.  Matthews 


SUMMARY 

An  analysis  is  made  of  the  effects  of  compressibility  on  the 
pressure  coefficients  about  several  bodies  of  revolution  by  com¬ 
paring  experimentally  determined  pressure  coefficients  with 
corresponding  pressure  coefficients  calculated  by  the  use  of  the 
linearized  equations  of  compressible  flow.  The  results  show 
that  the  theoretical  methods  predict  the  subsonic  pressure - 
coefficient  changes  over  the  central  part  of  the  body  but  do  not 
predict  the  pressure-coefficient  changes  near  the  nose.  Ex¬ 
trapolation  of  the  linearized  subsonic  theory  into  the  mixed, 
subsonic-supersonic  flow  region  fails  to  predict  a  rearward 
movement  of  the  negative  pressure-coefficient  peak  which  occurs 
after  the  critical  stream  Mach  number  has  been  attained.  Two 
equations  developed  from  a  consideration  of  the  subsonic  com¬ 
pressible  flow  about  a  prolate  spheroid  are  shown  to  predict, 
approximately,  the  change  with  Mach  number  of  the  subsonic 
pressure  coefficients  for  regular  bodies  of  revolution  of  fineness 
ratio  6  or  greater. 

INTRODUCTION 

A  number  of  papers  have  been  published  concerning  the 
theoretical  aspect  of  the  effects  of  compressibility  on  the 
flow  over  bodies  of  revolution  (refs.  1  to  4).  In  the  present 
investigation  these  theoretical  methods  are  applied  to  the 
analysis  of  experimental  data.  Such  an  analysis  should 
contribute  to  the  basic  knowledge  of  subsonic  three- 
dimensional  flow. 

Two  prolate  spheroids  of  fineness  ratios  6  and  10,  an  ogival 
body,  and  a  prolate  spheroid  with  an  annular  bump  near  the 
nose  were  tested  in  this  investigation.  The  experimental 
pressures  about  the  two  prolate  spheroids  are  compared  with 
the  pressures  computed  by  the  linearized  compressible-flow 
theory.  Several  relations  developed  from  theoretical  con¬ 
siderations  of  the  flow  about  a  prolate  spheroid  are  presented 
for  correcting  the  incompressible  pressure  coefficients  of 
regular  bodies  of  fineness  ratios  6  to  10  for  the  effects  of 
compressibility  in  the  subcritical  flow  range.  Results 
obtained  from  these  relations  are  also  compared  with  corre¬ 
sponding  experimental  pressure  coefficients. 

8YMBOLS 

6  maximum  radius  of  body 

Cm  normal-force  coefficient  based  on  plan-form  area  of 
ellipse 


f  fineness  ratio  of  body,  l/2b 

l  total  length  of  body  (see  fig.  1 ) 

Mc,  critical  Mach  number 

M0  free-stream  Mach  number 

Pi  local  static  pressure 

Po  free-stream  static  pressure 


pressure  coefficient, 


Po 


pV> 

£d 

r  local  radius  of  body 

S(x)  cross-sectional  area  of  body  of  revolution 

component  of  local  velocity  parallel  to  free  stream 


U 


V 

w 

x 


free-stream  velocity 
component  of  local  velocity  in  vertical  plane  perpen¬ 
dicular  to  free  stream 
total  local  velocity 

component  of  local  velocity  perpendicular  to  u  and  o 
coordinate  along  major  axis  of  body 
angle  of  attack 


y  ratio  of  specific  heat  at  constant  pressure  to  specific 

heat  at  constant  volume 
p  density 

<P  velocity  potential 

ft,  f,  w  ellipsoidal  coordinates  (see  ref.  5) 

Subscripts: 

c  compressible  value 

i  incompressible  value 

cr  critical  value 

st  incompressible  value  of  flow  about  hypothetical 

stretched  body 

MODELS 

Sketches  of  the  bodies  of  revolution  tested,  which  show 
the  locations  of  the  pressure  orifices  and  other  pertinent 
details,  are  presented  in  figure  1.  The  ordinates  of  the 
typical  transonic  or  ogival  body  and  the  prolate  spheroid 
with  an  annular  bump  are  given  in  table  I.  The  ordinates 
of  the  section  of  the  sting  support,  which  is  a  part  of  the 
body  of  revolution,  are  those  of  a  prolate  spheroid  of  fineness 
ratio  6.  The  same  support  was  used  for  each  body.  The 
couplings  used  to  change  the  angle  of  attack  were  mounted 
in  the  sting  11  inches  downstream  from  the  end  of  the  body. 
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TABLE  I.— ORDINATES  OF  THE  TYPICAL  TRANSONIC 
BODY  AND  OF  THE  ANNULAR  BUMP  PROLATE  8PHEROID 


Transonic  body 

Prolate  spheroid  with 
annular  bump 

rfl,  percent 

rll,  percent 

rll,  percent 

rll,  percent 

0.00 

0.000 

0.00 

0.000 

.60 

.  462 

.75 

1.437 

.76 

.686 

1.25 

1.854 

1.26 

-KM 

2.50 

2.604 

2.60 

1.445 

5.00 

3.631 

5.  no 

2.409 

10.00 

5.000 

10.00 

3.940 

12.50 

5.560 

20.00 

6.180 

12.91 

5.684 

30.00 

7.  480 

13.33 

5.873 

40.00 

8. 121 

14.  16 

6.408 

60.00 

8.333 

15.20 

7.230 

60.00 

8.162 

16.24 

7.800 

70.00 

7.636 

17. 

8.  190 

75.00 

7. 215 

17.91 

8.230 

18. 74 

8.250 

20.00 

8.230 

23. 32 

8.220 

26.00 

8.000 

27.08 

7.810 

28.31 

7.663 

28. 74 

7.620 

29.  16 

7.605 

30.00 

7.640 

35.00 

7. 952 

50.00 

8.290 

55.00 

8.333 

65.00 

8.290 

75.00 

7.  952 

j-bump- 

— i — — — — -Le“ — d - 1 - — 

*  c 
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(a)  Prolate  spheroid  ;/=  6. 

(b)  Prolate  spheroid;  /— 10. 

(c)  Typical  transonic  body. 

(d)  Prolate  spheroid  with  annular  bump. 

(e)  Angular  locations  of  orifices. 

Figure  I. — Profiles  of  bodies  tested. 

Except  at  the  first  three  stations  indicated  in  figures  1  (a) 
to  1  (d),  the  pressure  orifices  were  located  around  the  body  as 
shown  in  figure  1  (e).  These  orifices  were  spaced  15°  apart 
on  one  side  of  the  body  in  order  to  obtain  a  fairly  accurate 
normal-force  coefficient  upon  integration  of  the  pressure 
coefficients.  The  orifice  at  the  first  station  was  located  in 
the  nose.  The  orifices  at  the  next  two  stations  were  located 
at  90°  intervals  around  the  body.  The  pressure  orifice 
openings  were  0.010  inch  in  diameter. 

TESTS 

The  pressures  about  the  bodies  were  measured  in  the 
Langley  8-foot  high-speed  tunnel  through  the  Mach  number 
range  0.3  to  0.95.  The  angle-of-attack  ranges  were  0°  to 
7.7°  for  the  regular  bodies  and  0°  to  2°  for  the  prolate 
spheroid  with  an  annular  bump.  The  pressures  were  re¬ 
corded  by  photographing  a  10-foot  100-tube  manometer 
board  filled  with  acetylene  tetrabromide. 


The  free-stream  pressures  and  Mach  numbers  were  deter¬ 
mined  from  an  empty-tunnel  calibration  based  on  the 
pressures  at  an  orifice  located  4  feet  upstream  of  the  model. 

Several  preliminary  plots  of  local  pressure  coefficients  as 
functions  of  free-stream  Mach  number  showed  considerable 
scatter  for  Mach  numbers  less  than  0.5,  probably  because 
of  the  difficulty  of  reading  the  small  pressure  differences  and 
because  of  the  possibility  that  the  tunnel  was  not  held  at 
each  Mach  number  a  sufficient  length  of  time  to  insure  com¬ 
plete  settling  of  the  manometer  liquid.  Because  of  this 
scatter,  it  was  necessary  to  neglect  the  pressure  coefficients 
below  A/o=0.5  in  extrapolating  the  pressure-coefficient  curves 
to  a  stream  Maeh  number  of  zero.  The  data  used  in  the 
analysis  in  this  investigation  were  picket!  from  the  extra¬ 
polated  curves. 

For  the  tests  reported  herein,  the  Reynolds  number  varies 
from  approximately  2,700,000  per  foot  at  A/0=0.40  to 
3,950,000  at  A/0=0.94. 

The  wall  interference  may  be  approximately  determined 
by  using  the  equations  of  reference  6.  Since  the  corrections 
were  small,  they  were  not  applied  to  the  pressures  in  the 
figures  which  present  experimental  data  alone;  however,  the 
corrections,  even  though  small,  were  applied  to  the  experi¬ 
mental  data  used  for  the  comparisons  between  the  theoretical 
and  the  experimental  values. 

THEORETICAL  METHODS 

The  theoretical  subsonic  pressures  about  a  prolate  spheroid 
may  be  computed  by  applying  the  Prandtl-Glauert  correction 
to  the  incompressible  potential-flow  equations  in  the  manner 
suggested  in  reference  7.  In  this  solution  of  the  linearized 
form  of  the  equations  for  compressible  flow,  the  body  is 
stretched  in  the  free-stream  direction  by  the  factor  1/0;  the 
induced  velocity  components  u—U,  v,  and  w  about  the 
stretched  body  are  computed  by  potential-flow  methods 
(for  prolate  spheroids,  see  ref.  5) ;  and  the  induced  velocities 
u—  U,  v,  and  w  are  corrected  by  the  factors  1/0J,  1/0,  and  1/0, 
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respectively.  The  corrected  velocities  are  the  compressible 
velocities  at  the  corresponding  points  on  the  original  body. 
The  following  formula,  as  is  shown  in  appendix  A,  is  the 
result  of  the  application  of  this  method  to  the  flow  over 
prolate  spheroids: 


sin*  w  sin2  a„ 


+(#»„  sin2  o>  sin2  <*•>  +  — 

~(~G7,)  ~2Kb“  8in2  “  9in2  «"]} 


(1) 


where 


F, i=  y  1  —  m2  cos  a„  —  m\  1  —  e,,2  cos  u  sin  a,, 
G„—l—e„2ii2 

H ,<=  >  1  —n2Katl  cos  a„— —*u2Kbst  cos  co  sin  a,, 


2e„ 


1  -eJ 

l-e.,* 

lo„  1  +<’« 

\  2#,? 

,0gl-e„ 

)  1-e,,1 

l+«„\ 

2e„(l~2el 

l-e.t) 

1 

f.t 


-V'-? 


tan  a„  =  0  tan  a 


with  a  function  of  log  0  and/  rather  than  with  1/0  as  in  two- 
dimensional  flow.  Equation  (3)  may  be  reduced  to  the  form 

—  =  1  -4-—^  which  is  presented  in  reference  8. 

P,  1-log  2/ 

Another  and  easier  method  of  obtaining  an  approximate 
solution  of  the  linearized  equations  for  very  thin  bodies  may 
be  found  in  references  2  to  4.  This  method  consists  of  in¬ 
tegrating  an  approximated  source-sink  distribution  to  ob¬ 
tain  the  induced-velocity  ratios  from  which  the  pressure 
coefficients  may  be  computed.  Since  the  source-sink  dis¬ 
tribution  is  approximated  by  the  derivative  of  the  cross- 
sectional  area  with  respect  to  the  length  of  the  body,  this 
method  is  more  generally  applicable  to  bodies  of  revolution 
than  is  the  method  of  applying  the  Prandtl-Glauert  correc¬ 
tion  to  the  exact  incompressible-flow  solution.  It  is  shown 
in  appendix  A  that,  for  prolate  spheroids  at  zero  angle  of 
attack,  this  method  gives  the  following  result : 


P  =— 
'  f 


h<-;> 

-W( 


+  02 


-  +  02- 

l2+P  l 2 


log 


-f+Vr+'F 


(4) 


Two  approximate  forms  which  show  the  effects  of  com¬ 
pressibility  can  be  obtained  from  equation  (4)  by  considering 
(a)  the  difference,  ( Pe—Pt ),  and  (b)  the  ratio  PJPt  of  the 
compressible  and  the  incompressible  values.  These  two 

0V2 

P 


relations  may  be  reduced  to  the  following  forms  when 


is  considered  small  with  respect  to 


(‘-O' 


or 


I2' 


The  pressure  coefficients  may  be  computed  from  the  fol¬ 
lowing  relation: 


Because  of  the  nature  of  the  transformation,  equation  (1) 
does  not  hold  for  large  angles  of  attack  (that  is,  where 
a=sin  a  ceases  to  be  a  fair  approximation)  or  for  bodies 
of  small  fineness  ratio. 

The  compressibility  effects  indicated  by  application  of  the 
linearized  theory  of  compressible  flow  to  prolate  spheroids 
are  not  apparent  from  equations  (1)  and  (2).  The  effects 
may  be  shown  simply  for  the  special  case  of  the  center  of  a 
prolate  spheroid  at  zero  angle  of  attack.  As  shown  in 
appendix  B,  the  following  relation  is  obtained: 

p-<=(i+ r  p- log  2  j  -] 

Pi  v  + 1  -log  2f)  L/2-0!Oog  2/- log  0)J  (3) 

Thus,  the  theoretical  solution  indicates  that  the  ratio  of  the 
compressible  pressure  coefficient  to  the  incompressible  pres¬ 
sure  coefficient  on  bodies  of  revolution  will  vary  conformably 


P-Pi 

2  log  0 

r 

(5) 

MX 

II 

+ 

log0 

1-log  2/ 

(6) 

Both  relations  indicate  that  the  effect  of  compressibility  on 
the  subsonic  flow  about  a  body  of  revolution  at  any  given 
Mach  number  is  to  lower  the  pressure  coefficients  over  a 
large  part  of  the  body.  These  relations  for  the  effect  of 
compressibility  are  in  accord  with  similar  equations  pre¬ 
sented  in  references  2  and  3. 

RESULTS  AND  ANALYSIS 

COMPARISON  OP  EXPERIMENTAL  AND  THEORETICAL  PRESSURE 
DISTRIBUTIONS 

The  local  pressure-coefficient  distributions  are  presented 
in  figures  2  to  6  for  various  values  of  free-stream  Mach 
number.  Figures  7  to  9  are  replots  of  some  of  the  data  of  the 
preceding  figures  corrected  for  wall  interference,  together 
with  results  of  the  theoretical  calculations  by  means  of 
equations  (2)  and  (4).  Figures  2  to  6  show  a  decrease  in  the 
experimental  pressures  over  the  central  part  of  the  body 
with  increasing  Mach  number,  as  predicted  by  equation  (5). 
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However,  figures  7  to  9  indicate  that  the  linearized  theory  effect  of  compressibility  on  the  experimental  pressure  oo- 

predicts  a  decrease  in  the  pressures  over  the  entire  body,  efficients  is  approximately  to  rotate  the  pressure-coefficient 

whereas  the  experimental  data  show  that  a  point  on  the  body  distributions  about  the  point  at  which  the  incompressible 

exists  ahead  of  which  the  pressures  increase  rather  than  pressure  coefficient  is  zero.  The  actual  point  about  which 

decrease.  (See  also  figs.  2  to  6.)  The  lack  of  agreement  of  the  rotation  may  be  considered  to  take  place  shifts  its 

the  linearized  theory  with  the  experimental  results  near  the  location  from  slightly  downstream  of  the  stream-pressure 

nose  of  the  body  is  to  be  expected  because  of  the  assumptions  point  on  the  top  of  the  body  to  slightly  upstream  of  the 

made  in  its  derivation.  It  might  be  pointed  out  that  the  stream-pressure  point  on  the  bottom  of  the  body. 


(a)  Prolate  spheroid; /=6;  Af„=0.916.  (b)  Prolate  spheroid;  /=  10;  A/«,=0.947. 

(e)  Typical  transonic  body;  0.898.  (d)  Prolate  spheroid  with  annular  bump;  Af ,,=0.634. 

Proms  3. — Experimental  pressure  distributions  over  several  bodies  of  revolution  at  sero  angle  of  attack. 


Pressure  coefficient,  P  Pressure  coefficient 
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(a)  «=0°. 

Fioubk  3. — Experimental  pressure  distributions  over  a  prolate  spheroid  of  fineness  ratio  6  at  several  angles  of  attack. 
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(b)  «=90°. 

Figure  3. — Continued. 
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(0)  a—  180°. 
Figure  3. — Concluded. 
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Pressure  coefficient,  P  Pressure  coefficient,  P 


REPORT  1155 — NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUTICS 


— —Ill— — 

BHGMMHH 


a --2.3"-  -Wcr  =  0942 


a  =37°;  Mcr  = 0.928 


a  *  5.9  ;  Afc-  =  0.93.5 


a  =  7.6  ;  Mgr 5  0.925 


50  0  I' 

Percent  distance  from  nose,  VI 


(a)  w=0°. 

Fioube  4. — Experimental  pressure  distributions  over  a  prolate  spheroid  of  fineness  ratio  10  at  several  angles  of  attack. 
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Figure  5. — Experimental  pressure  distribution  over  a  typical  transonic  body  at  several  angles  of  attack 
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As  the  flow  approaches  and  exceeds  the  critical  stream 
Mach  number,  a  further  change  in  the  pressure  distributions 
occurs.  This  change  of  shape  (figs.  2  to  5)  is  essentially  a 
rearward  movement  of  the  negative  pressure  peak.  The 
nature  of  this  change  is  emphasized  in  the  {dots  for  Af0= 0.950 
of  figures  7  and  8.  These  figures  show  that  the  linearized 
theory  does  not  predict  the  shift  in  peak  pressures  which 
occurs  as  the  flow  becomes  supercritical. 

The  rearward  shift  of  negative  pressure  peaks  which  occurs 
on  the  top  of  the  body  (figs.  3  (a),  4  (a),  and  5  (a))  seems 
to  be  changed  to  a  forward  shift  on  the  bottom  of  the  body 
(figs.  3  (c),  4  (c),  and  5  (c)).  It  is  reasonable  to  assume 
that  part  or  all  of  this  forward  movement  of  the  bottom 
negative  pressure  peak  may  be  explained  by  the  positive 
pressure  field  which  exists  ahead  of  the  under  part  of  the 
sting  support. 

A  comparison  of  figures  7  and  9  shows  that  the  linearized 
theory  gives  better  results  for  the  body  of  larger  fineness 
ratio.  The  pressures  about  the  prolate  spheroid  of  fineness 
ratio  10  are  in  better  agreement  with  theory  even  for  the 
stream  Mach  number  of  0.950  than  are  the  pressures  about 


the  body  of  fineness  ratio  6.  It  may  also  be  observed  that 
the  theoretical  pressures  about  the  prolate  spheroid  of  fine¬ 
ness  ratio  10,  which  are  calculated  by  the  two  different 
methods,  are  in  excellent  agreement;  thus,  these  results 
show  that,  for  bodies  of  fineness  ratios  of  10  or  greater,  the 
simpler  method  of  computing  pressures  presented  in  refer¬ 
ences  2  to  4  is  fairly  reliable. 

INFLUENCE  OF  CHANGING  NOSE  SHAPE 

The  effects  of  changing  the  shape  of  the  nose  of  a  body 
are  seen  by  comparing  figures  2  (a)  and  3  with  figures  2  (c) 
and  5.  The  incompressible  pressure  distribution  is  changed 
as  may  be  expected.  However,  the  nature  of  the  effect  of 
compressibility  is  the  same  for  this  body  as  for  the  prolate 
spheroid  of  fineness  ratio  6.  The  incremental  pressure 
changes  are  almost  the  same,  and  the  rotation  and  shifts  of 
pressure  peaks  are  very  similar  for  both  bodies.  This  com¬ 
parison  shows  that  the  effects  of  compressibility  do  not 
depend  to  a  great  extent  on  body  shape  so  long  as  the  body 
does  not  depart  from  the  specifications  required  for  the 
application  of  the  linearized  equations. 


Fiouee  7. — Comparison  of  experimental  and  theoretioal  pressures  over  a  prolate  spheroid  of  fineness  ratio  6  at  zero  angle  of  attack. 
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Figube  8. — Continued. 
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Figure  9. — Comparison  of  experimental  and  theoretical  pressures  of  a  prolate  spheroid  of  fineness  ratio  10  at  zero  angle  of  attack. 


INFLUENCE  OP  FINENESS  RATIO 

The  influence  of  fineness  ratio  on  the  effects  of  compressi¬ 
bility  may  be  observed  by  comparing  figures  2  (a)  and  3  with 
figures  2  (b)  and  4.  These  figures  show  that  increasing  the 
fineness  ratio  reduces  the  changes  in  pressure  caused  by 
varying  the  stream  Mach  number.  This  effect  is  predicted 
by  the  linearized  theory  in  equation  (5).  It  may  also  be 
observed  that  the  pressure  peaks  are  less  prominent  and  do 
not  shift  their  location  to  the  extent  found  for  the  bodies  of 
lower  fineness  ratio.  The  changes  in  the  shape  of  the  pres¬ 
sure  distributions  are  also  reduced  and  comparable  changes 
occur  at  higher  Mach  numbers.  The  delay  in  the  change  of 
the  shape  of  the  pressure  distribution  is  demonstrated  by 
comparing  figures  7  and  0  at  M0=0.95.  For  the  prolate 
spheroid  of  fineness  ratio  6  a  marked  change  in  the  pressure 


distribution  has  already  occurred,  whereas  for  the  body  of 
fineness  ratio  10  the  shape  of  the  pressure-distribution  curve 
is  almost  the  same  as  at  lower  Mach  numbers.  A  considera¬ 
tion  of  the  observed  effects  of  increasing  the  fineness  ratio 
indicates  that  such  a  change  definitely  reduces  the  effects 
of  compressibility. 

INFLUENCE  OF  ANGLE  OF  ATTACK 

It  may  be  shown  by  the  use  of  the  linearized  theory  that, 
at  least  to  a  first  approximation,  the  lift  and  moment  forces 
on  a  body  of  revolution  are  not  affected  by  changes  in  Mach 
number.  (See  ref.  4.)  The  validity  of  this  prediction  is 
demonstrated  in  figure  10  which  shows  that  the  variation  of 
the  normal-force  coefficient  with  Mach  number  is  small  for 
both  the/=10  and/=6  prolate  spheroids. 


Figure  10. — Theeffect  of  compressibility  on  the  normal-force  coefficients 
over  the  forward  half  of  two  prolate  spheroids  of  fineness  ratios  6  and  10. 

INFLUENCE  OF  AN  ANNULAR  BUMP 

A  study  of  the  effects  of  compressibility  on  the  velocities 
about  an  infinitely  long  body  containing  surface  waves  (ref. 
9)  shows  that  these  effects  become  two-dimensional  in  nature 
when  the  length  of  the  surface  waves  becomes  small  with 
respect  to  the  body  radius.  Since  an  annular  bump  on  a 
body  of  revolution  approximates  these  conditions,  the  flow  over 
such  a  bump  may  also  be  expected  to  show  two-dimensional 
effects.  An  examination  of  figures  2  (d)  and  6  shows 
that  the  range  of  pressure  coefficients  found  in  the  flow  over 
a  prolate  spheroid  with  an  annular  bump  is  of  the  same  order 
as  that  found  in  two-dimensional  flow.  The  two-dimensional 
nature  of  the  flow  over  an  annular  bump  is  further  demon¬ 
strated  by  comparing  the  pressure  coefficients  with  the  Von 
K&rm6n  relationship  (ref.  10)  for  the  effects  of  compressibility 
on  two-dimensional  flow  (fig.  11).  Figure  11  shows  fair 
agreement  between  the  Von  K£rm£n  relation  and  the  experi¬ 
mental  relationships  for  those  regions  of  the  body  where  the 
flow  does  not  separate  and  the  slope  of  the  body  is  reasonably 
small;  namely,  the  8.33-,  11. 6-,  13.6-,  10.5-,  17.9-,  and  19.8- 


percent  stations.  The  15-percent  station  is  highly  irregular 
and  cannot  be  explained  by  either  two-  or  three-dimensional 
theories.  The  other  stations  are  severely  affected  by  separa¬ 
tion  phenomena.  The  Von  Kdrmdn  relation,  however,  fails 
to  explain  the  phenomena  once  the  critical  speed  is  exceeded. 

CORRECTION  OF  INCOMPRESSIBLE  PRESSURE  DISTRIBUTIONS  FOB  THE 
EFFECTS  OF  COMPRESSIBILITY 

Equations  (5)  and  (6)  suggest  that  an  incompressible  pres¬ 
sure  distribution  might  be  corrected  for  the  effects  of  com¬ 
pressibility  by  considering  a  pressure-increment  type  of 
function  such  as  Pc— Pt  or  a  rate-of-increase  type  of  function 
such  as  PJPt .  In  order  to  show  whether  the  effects  of  com¬ 
pressibility  may  be  expressed  by  such  functions,  a  number 
of  the  pressures  over  the  regular  bodies  at  zero  angle  of 
attack  have  been  plotted  in  figure  12  in  terms  of  PJPt  and 
Pc— Pt  against  xfl  and  M0.  Tunnel-wall  corrections  have 
been  omitted,  but  the  omission  does  not  affect  the  conclu¬ 
sions.  An  examination  of  both  functions  shows  that,  except 
at  supercritical  Mach  numbers,  the  values  of  PJPi  and 
Pe—P,  are  roughly  constant  between  the  25-percent  and  the 
50- percent  stations.  Over  the  forward  part  of  the  body,  the 
values  are  more  variable. 
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The  PJP,  function  becomes  discontinuous  in  the  neigh¬ 
borhood  of  Pi— 0.  This  behavior  may  be  attributed  to  the 
fact  that  the  pressure  coefficient  is  zero  at  the  incompressible 
stream-pressure  point  and,  since  one  of  the  effects  of  com¬ 
pressibility  is  to  shift  the  stream-pressure  point,  discon¬ 
tinuities  may  be  expected  in  the  neighborhood  of  this  point. 
However,  since  the  pressures  in  this  region  are  small,  a  wide 
variation  in  Pe/Pt  may  be  permissible  without  serious  error 
in  the  corrected  results. 

The  Pc—P,  correction  may  also  be  expected  to  become 
irregular  in  the  region  of  the  nose.  The  experimental  curves 
show  that  this  function  changes  sign  in  the  neighborhood  of 
the  stream-pressure  point  so  that  any  correction  function  of 
this  type  should  Liclude  the  position  on  the  body.  However, 
such  a  function  cannot  be  obtained  from  the  linearized 
method  as  this  method  does  not  indicate  the  change  of  sign 
shown  in  the  experimental  data. 

The  experimental  values  of  Pe—Pt  and  Pf/P<  at  the  cen¬ 
ters  of  the  regular  bodies  are  compared  with  equations  (5) 
and  (6)  in  figure  13  in  order  to  show  the  validity  of  the  pre¬ 
diction  of  the  effect  of  compressibility  by  the  linearized 
potent’'!-!.  >w  theory.  It  is  observed  that  equation  (6) 


Figure  13. — Theoretical  compressibility  correction  functions  com¬ 
pared  with  experimental  results.  a=0°;  -j=S0  percent. 


within  its  limitations  predicts  the  effects  of  compressibility 

Pc  1 

for  three-dimensional  flow  whereas  the  relation  pC==^' 


which  is  used  to  predict  the  compressibility  effects  of  two- 
dimensional  flow,  does  not.  It  may  also  be  observed  that 
equations  (5)  and  (6)  predict  the  effects  of  compressibility 
with  about  the  same  degree  of  accuracy. 

The  correction  functions  are  applied  to  several  incompres¬ 
sible  pressure-coefficient  distributions  in  figure  14,  which 
are  compared  with  the  corresponding  experimental  distribu¬ 
tions.  It  is  shown  in  figure  14  (a)  that  increasing  the  fineness 
ratio  of  the  prolate  spheroid  from  6  to  10  or  reducing  the 
bluntness  of  the  nose,  which  is  the  essential  difference  between 
the  ogival  body  and  the  prolate  spheroid,  extends  the  region 
of  the  body  for  which  corrections  can  be  made  from  the 
20-percent  station  for  the  prolate  spheroid  of  fineness  ratio  6 
forward  at  least  to  the  10-percent  station  for  the  sharper-nose 
bodies.  The  PjPi  function  expresses  the  effect  of  compres¬ 
sibility  more  accurately  in  the  vicinity  of  the  nose  than  does 
the  Pc—P t  function.  This  result  is  to  be  expected  since  one 
of  the  effects  of  compressibility  already  noted  is  the  rotation 
of  the  pressure  distribution,  which  is  accounted  for  by  the 
PJPi  expression  but  not  by  the  Pc—Pt  expression. 

The  increasing  error  which  results  from  increasing  the 
stream  Mach  number  is  shown  in  figure  14  (b).  At 
M0=  0.800,  the  incompressible  pressure  coefficients  about  the 
fineness  ratio  6  prolate  spheroid  may  be  corrected  with  a 
fair  degree  of  accuracy  as  far  f<  u  ard  as  the  5-percent  station. 
As  the  Mach  number  ii  ° reuses,  the  divergence  between  the 
corrected  values  and  the  experimental  values  in  the  region 
of  the  nose  increases  and,  with  still  greater  Mach  numbers, 
tends  to  spread  toward  the  center.  At  M0—  0.940,  which  is 
supercritical  for  the  prolate  spheroid  of  fineness  ratio  6,  the 
correction  formulas  are  still  applicable  at  the  center,  so  that 
successful  extrapolation  of  the  linearized  theory  into  the 
supercritical  range  is  found  to  depend  on  the  section  of  the 
body  to  which  the  extrapolation  is  applied. 

As  may  be  expected,  the  success  of  the  linearized  theory  in 
expressing  the  effects  of  compressibility  decreases  as  the 
angle  of  attack  increases.  The  principal  reason  for  this 
result  is  that  an  angle  of  attack  involves  a  pressure  peak  on 
the  forepart  of  the  top  of  the  body,  which  moves  rearward 
when  the  stream  Mach  number  approaches  and  exceeds  the 
critical  value  for  the  body.  Since  the  correction  formulas 
either  rotate  or  translate  the  incompressible  pressure  distri¬ 
bution,  they  cannot  express  this  change  in  the  shape  of  the 
pressure  distribution.  This  phenomenon  is  demonstrated 
in  figure  14  (c),  which  presents  a  comparison  of  the  corrected 
pressure-coefficient  distributions  and  the  experimental  dis¬ 
tributions  of  tbe  flow  about  the  prolate  spheroid  of  fineness 
ratio  6  at  several  angles  of  attack.  Even  though  the  shift  of 
the  peak  pressure  is  not  accounted  for  in  the  correction  for¬ 
mula,  the  corrected  distributions  are  not  seriously  in  error 
at  the  peaks  and  the  agreement  improves  over  the  midportion 
of  the  body.  Thus,  if  some  error  is  permissible,  these  for¬ 
mulas  may  be  applied  for  angles  of  attack  as  high  as  7°  or  8°. 
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(c)  /=  6;  prolate  spheroid;  A/o=0.900. 
Figure  14. — Concluded. 


Figure  14  (c)  indicates  that  equation  (6)  does  not  correct 
as  satisfactorily  over  the  central  parts  of  the  body  at  angles 
of  attack  as  equation  (5).  This  lack  of  agreement  is  due  to 
the  compressibility  effect  on  the  lift  forces.  It  has  already 
been  shown  that  the  lift  forces  are  not  much  affected  by 
compressibility;  hence,  the  increments  of  the  pressure 
coefficients  due  to  compressibility  are  about  the  same  for 
the  top  and  bottom  of  the  body.  Since  the  absolute  values 
of  the  pressure  coefficients  are  less  on  the  bottom  of  the 
body  and  greater  on  the  top  than  if  the  lift  forces  had  not 
been  present,  equation  (6)  will  overcorrect  the  pressure 
coefficients  on  the  top  and  undercorrect  those  on  the  bottom. 
The  same  reasoning  shows  that  equation  (5),  which  gives 
a  constant  increment  over  the  entire  body,  will  express  the 
compressibility  effect  with  an  angle  of  attack  better  over 
the  central  part  of  the  body  to  which  it  applies  than  will 
equation  (6). 

CONCLUSIONS 

The  results  of  the  tests  made  on  several  bodies  of  revolu¬ 
tion  have  shown  the  following  effects  of  compressibility  on 
three-dimensional  flow: 

1.  In  general,  the  compressibility  effect  is  to  increase  the 
pressure  differences  over  a  body  of  revolution.  The  pressure 
distributions  are  approximately  rotated  about  a  point  near 
stream  pressure  and  the  negative-pressure  peaks  are  moved 
rearward. 

2.  The  linearized  modification  of  the  compressible 
potential-flow  equation  will  predict  the  pressures  over  the 
central  part  of  the  body  but  will  not  predict  the  changes  in 
pressure  ahead  of  the  stream-pressure  point  nor  will  it  pre¬ 
dict  the  change  in  shape  which  occurs  with  supercritical  flow. 


3.  The  correction  formu! 

Zj_i  ■  ft 

Pt  '-log  2/ 

and 


(where  Pe  and  P,  are  the  pressure  coefficients  for  compress¬ 
ible  and  incompressible  flow,  respectively,  /  is  the  fineness 
ratio,  and  /3 = yl  —~MJ  in  which  M0  is  the  Mach  number) 
may  be  used  approximately  to  correct  incompressible-flow 
pressures  over  the  central  part  of  streamline  thin  bodies  of 
revolution;  the  errors  will  increase  as  the  supercritical  Mach 
number  is  approached  and  exceeded.  Since  PJPi  rotates 
the  pressure  distribution,  this  correction  is  better  to  use  at 
zero  angle  of  attack;  however,  the  form  Pc— P<  expresses 
the  effects  of  angle  of  attack  more  correctly  and  should  be 
used  when  an  angle  of  attack  other  than  zero  is  involved. 

4.  The  effects  of  compressibility  are  approximately  the  same 
for  various  bodies  of  the  same  fineness  ratio,  provided  the 
body  shape  satisfies  the  requirements  of  the  linearized  theory. 

5.  Increasing  the  fineness  ratio  tends  to  reduce  the  effects 
of  compressibility. 

6.  The  effects  of  compressibility  on  an  annular  protuber¬ 
ance  of  short  chord  on  a  body  of  revolution  tend  to  follow 
more  nearly  two-dimensional  laws  than  three-dimensional 
laws. 

7.  Lift  forces  and  moments  over  the  forward  part  of  the 
body  are  relatively  unaffected  by  compressibility. 

Langley  Aeronautical  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Langley  Field,  Va.,  November  6,  1951. 


APPENDIX  A 


DERIVATION  OF  THE  EQUATIONS  FOR  THE  COMPRESSIBLE  PRESSURE  COEFFICIENTS  OF  THE  FLOW  ABOUT  A  PROLATE 

SPHEROID 

The  solution  of  the  linearized  compressible-flow  equation  velocity  components  obtained  by  differentiating  the  potential 
for  a  prolate  spheroid  requires  a  derivation  of  the  relation  equation  are 

for  the  incompressible  velocities  about  the  body.  The  in-  , _ _ 

compressible  velocities  about  a  prolate  spheroid  are  defined  ^ _  1  ~m2  ^  cog  a  MV£lzzM*KLlZff)  g* oos  u 8jn  a  1 

bv  the  ootential  eouations  (riven  in  reference  5.  These  £  l  —e  n2  l—en‘ 


by  the  potential  equations  given  in  reference  5.  These 
equations  may  be  combined  and  written 

<p—Un{  cos  a  +  A  (cos  <*)m^  f  log 

B  (sin  «)  Vr=7Vf^  (|  log  fzr]-ptr)  cos  “  (7) 


-**)K 

U~  I-eV 


cos  u cos  a 


„2(1  —e1) 

4-  — r  -f  Kb  cos2a>  sin  a-\-Kt  sin2a>  sin  a 

1  — e  u 


The  values  of  the  constants  A  and  B  which  satisfy  the  w* _ mV(1  — M*)(l  —  e2)  rs  ■ 

required  boundary  conditions  are  U  1— c2m2  « sin  a>cosa 

u  u*(\—  __  . 


'-Is _ I  w  is±i 

fo2-  1  2  10g  fo-  I 


,Io2-J _ 1,  £o±l 

f„(f„2-l)  2,Ogf0-l 

where  f0  is  the  value  of  the  coordinate  which  represents  the 
body.  It  may  be  shown  that  the  eccentricity  of  the  ellipse 

€  =  , JL  where  l  and  26  are  the  lengths  of  the  major 

and  minor  axes  of  the  prolate  spheroid.  Since  the  fineness 


+  ~r — 5— y  Kb  sin  w  cos  a)  sm  a — A 1,  sin  tocos  a)  sm  a 
1  — e  u.‘ 


where 


Ka=  1- 


log(l^c)_  2c 


log! 

( 

,  2c 

u- 

*7 

1  l_e2 

log! 

(S) 

1- 

2e(l-2c2j 

1  ^ 

ratio  /  is  equal  to  l/2b,  e--^l ~Ji‘  The  incompressible  antj  u*,  t?*)  an,i  w*  are  the  velocity  components  in  a  coordi- 

- - -  nate  system  alined  with  the  /-axis  of  the  body.  These 

velocities  are  transformed  to  the  u,  v,  and  w  components  by  the  following  equations: 

■V 

u  v*  ,  r*  . 
jj—jf  cos  a  +  Sill  a 

B  U*  .  .  V* 

jj—  —  jj  sin  a+jr  cos  a  >-  (10) 

w  w* 

U~V 

* 

With  the  preceding  transformation,  the  velocity  equations  (8)  become 
p=  (\1  ~  (t2  cos  a — mV  1  — e2  cos  w  sin  a)(\l  — M2  K,  cos  a—  tt\l—  e2  Kb  cos  to  sin  a)+Kb  sin2«  sin2a 

p  j  _  _  _  _ 

jj  —  ~  l  _gy»(V  1~M2  sin  a+fty  l  —e2  cos  a  cos  a)(V  1  — p2  Aa  cos  a~nfj\  —el  Kb  cos  te  sin  <*)-f  Kb  sin2ai  cos  a  sin  a  (j  j) 


^  sin  ( V 1  —  M2  A.  cos  a— mV1  ~«2  cos  u  sin  a)— Kb  sin  «  cos  «  sin  a 
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These  equations  may  be  rewritten  more  simply  by  setting 


F—  yrl  —  n*  cos  a — n->J\ — e1  cos  1a  sin  a 
Fi=  Vl  — M*  sin  a+in/l—e2  cos  u>  cos  a 
H—  Vl  —  m1  Kt  cos  a — mVI — e*  K*  cos  “  9*n  a 


<7=  1-eV 


—e2  sin  a 


u  FH .  „  .  .  .  » 

7j—  'q~  +K*  sra*»  sura 


v  FiH  ,  v  ■  , 

U~ - 6 — '  ■**  9m  “  cos  «  sm  a 


w  FiH 


w  l  2*a  rr  •  • 

jj~ - -q - Kb  sm  ui  cos  to  sm  a 


The  method  of  correcting  for  compressibility  discussed  in 
the  text  can  now  be  approximately  applied  by  increasing  the 
fineness  ratio  by  1//J  and  reducing  the  tangent  of  the  angle 
of  attack  by  the  factor  0.  Thus, 


tan  tan  a 


where  e„  and  a„  are  the  eccentricity  and  the  angle  of  attack 
of  the  stretched  body.  Although  this  stretched  body  differs 
slightly  from  the  properly  stretched  body,  the  approximation 
is  very  close  for  large  fineness  ratios  and  small  angles  of 
attack.  It  may  be  shown  that  the  properly  stretched  body 
is  an  ellipsoid  having  three  unequal  axes;  however,  under  the 
present  restrictions  the  two  minor  axes  are  very  nearly 
equal,  so  that  only  small  errors  will  be  caused  by  the  above 
approximation  of  the  stretched  body. 

The  induced  velocities  in  compressible  flow  are  now  de¬ 
termined  by  substituting  a„  and  e„  in  equations  (11)  and 


li  V  2V 

multiplying  the  resulting  velocity  increments  jj—  1 ,  and  jy 

by  the  factors  l/0a,  1//3,  and  1//3,  respectively,  or 


(u)ri(ul 


The  pressure  coefficients  for  compressible  flow  may  be  com¬ 
puted  from  the  velocities  by  the  following  formula: 


where 


(«)>(«>: 


Combination  of  equations  (12),  (15),  and  (16)  yields 


,  V2  1  (.  H.,2  v  i  ■  i  -2 

1  ~ip=  Ji  \  1  8in  "  mb*  a.,- 

(?-  o  [(■  -+^)'- 


(%?)2_2^*1  8in2  “ 8in2  “■•]} 


where  the  subscript  st  is  used  to  indicate  that  the  various 
functions  so  identified  are  based  on  the  values  of  e  and  a 
associated  with  the  stretched  body.  (See  eqs.  (13)  and  (14).) 
A  simpler  first  approximation  may  be  obtained  by  consid¬ 
ering  the  approximate  relation 


m. 


(18) 


A  simpler  equation  may  be  developed  for  the  pressures 
over  a  prolate  spheroid  at  zero  angle  of  attack  by  consider¬ 
ing  the  method  of  approximate  source-sink  distributions 
described  in  references  2  to  4.  In  these  references,  it  is 
shown  that 


p  - _ 9  (u-u\  __L  f 1  S'(t)(x~t)dt 

rc~  zy  v  Jr2rJo[(x_l)1+^v2 


where  t  is  a  coordinate  along  the  major  axis  of  the  body  and 
S'(t)  is  the  derivative  of  the  cross-sectional  area  of  the  body 
with  respect  to  t.  For  a  prolate  spheroid, 


o/A  _.2  4rb2t2 

i>(t) — rr  ^  p 


from  which 


p_2 b2  f‘V 
e~  I-  Jo  [(*- 


i  (1-y)  (x_<) dt 


t)2+P2r2]il2 


After  integration  and  collection  of  terms 


L2V(1-T )+f*1Ti 

,-f+VRRIl 


“7  +\li+P*Ti 


APPENDIX  B 


REDUCTION  OF  PRESSURE-COEFFICIENT  FORMULAS  TO  OBTAIN  SIMPLE  FUNCTIONS  FOR  CORRECTING 
INCOMPRESSIBLE  PRES8URE  DISTRIBUTIONS  FOR  THE  EFFECTS  OF  COMPRESSIBILITY 


Two  functions  which  may  be  used  to  express  the  relation 
between  the  pressure  coefficients  in  compressible  and  incom¬ 
pressible  flow  are  the  ratio  and  the  increment  between  the 
two  coefficients  P„  and  Pc,  that  is,  PJPt  and  Pc—Pt.  Both 
functions  may  be  expressed  in  simple  equations  by  substitut¬ 
ing  the  pressure-coefficient  functions  for  the  midpoint  of  the 
body  into  both  the  ratio  function  and  the  increment  func¬ 
tion.  In  order  to  simplify  equation  (17)  let  n—0,  sin  a,,=/3a, 
and  K,t=  1 — k„  or  k,,=  1 —Katt.  Then, 

V2 

For  small  values  of  1—jp 

V 2 

P  —  i _ L_ 

1  U2 

Also,  for  large  values  of/ 

K^r*  2 

Hence 

Since  /3*=1  at  Af0=0 

Pe  1  [2fc.,-^-4aJ(l-/32)(/3  V-*„)] 

frj'  W^k? 

Since  4a*(l— ^a)OV— it,,)  is  small  compared  with  2 k,„  the 
term  containing  a  may  be  neglected ;  thus 


9  *«< 

Pc _ 1  k„  p2 

Pt  /S*  ki  2-kt 


(22) 


In  order  to  reduce  this  equation  to  previously  published 
forms  (refs.  4  and  8),  it  is  necessary  to  reduce  k: 


S2  B2 

Substituting  e„2= 1 — and  the  approximate  form  e„=  1 — ~p 

in  this  equation  gives 


.  3*  (log  2/— log  B~  1) 
'B2  log  2f-p2  log  B-f2 


l _ i°g  2/  1 

*l— log  2 f  —  f2  {  > 

These  equations  show  that  both  /t.,/01  and  kf  are  of  order  of 
magnitude  Ilf2  and,  therefore,  small  with  respect  to  2. 
Hence,  the  approximation  (see  eq.  (22)) 

pre2  ^  (25) 

is  valid. 

If  equations  (23)  and  (24)  are  used  in  equation  (25),  the 
following  equation  is  obtained: 

Ps=(1+  r _ /2~Iog  2/ _ 1  (26) 

Pi  v  + 1  -log  2fJ  Lf-B2  (log  2/— log  B)  J  K  > 

or  for  large  fineness  ratios 

P4“1  +  l-log2/  (27) 

which  may  be  changed  to  its  equivalent  form 

Pc— Pi— ^—yf-  (28) 


Equation  (20)  obtained  by  the  source-sink-distribution 
method  will  also  reduce  to  equations  (27)  and  (28)  for  the 
central  part  of  the  body. 
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